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Abstract 

The effect of strong magnetic field on the stability and gross properties of bulk 
as well as quasi-bulk quark matter has been investigated using the conventional 
MIT bag model. Both the Landau diamagnetism and the paramagnetism of quark 
matter has been studied. How the quark hadron phase transition is affected by the 
presence of strong magnetic field has also been investigated. The equation of state 
of strange quark matter changes significantly in strong magnetic field. It is also 
shown that the thermal nucleation of quark bubbles in compact metastable state 
of neutron matter is completely forbidden in presence of strong magnetic field. 
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1. Introduction 

It is generally believed that if the density of neutron matter at the core of 
a neutron star exceeds a few times normal nuclear matter density, a deconfining 
transition to quark matter may take place. Whatever be the order of this transi- 
tion, will convert a normal neutron star to a hybrid star with an infinite cluster 
of quark matter at the core and a crust of neutron matter. If the speculation 
of Witten [1] that a flavor symmetric strange quark matter (SQM) is the true 
ground state of matter at zero temperature and pressure is correct, then there 
is a possibility that the whole neutron star will be converted to a star of SQM, 
in which strange quarks are produced through the weak processes [2] and ulti- 
mately giving rise to a dynamical chemical equilibrium among the constituents. 
The stability of such bulk and quasi-bulk objects have been thoroughly studied 
using MIT bag model [3-6] as well as other confinement models, e.g. we have 
studied using Dynamical Density Dependent Quark Mass [2,7-9] (D3QM) model 
of confinement, both at zero and finite temperatures. In this context we should 
mention that the results of MIT bag model does not agree with the recent lattice 
Monte Carlo results. The latter predicts a weak first order or a second order QCD 
phase transition, whereas according to MIT bag model it is a strong first oredr 
transition. On the other hand, we have seen that the predictions of D3QM model 
more or less agree with the recent lattice Monte Carlo results. Since we do not 
know how to incorporate the effect of strong magnetic field in lattice gauge calcu- 
lations, we shall concentrate only on the MIT bag and D3QM models. To see the 
effect of strong magnetic field on second order quark-hadron phase transition we 
shall assume a metal-insulator type of transition in presence of strong magnetic 
field. The stability and gross properties of hybrid and strange stars have also been 
investigated [10-12]. Although the gross properties of such self-gravitating objects 
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do not differ significantly from a pure neutron star (except the M-R diagram), 
the presence of quark phase at the core may cause rapid cooling of the star and 
can also produce significant gamma ray bursts. The other important aspect is 
the superconducting nature of the quark matter core. Since quarks carry electric 
charges, unlike neutron matter which exhibits superfluidity in a neutron star, this 
particular phase as a whole behave like a superconducting fluid (and possibly type 
I; of course if there are very few protons present in neutron matter, they may 
form Cooper pairs and give rise to superconducting zones). As a result in such 
objects the magnetic field lines mainly pass through the crust region (if any), the 
consequence of which is the rapid Ohmic decay of the magnetic field. Therefore 
the surface magnetic field for such objects will be low enough compared to the 
corresponding pure neutron star. 

The other interesting part of such objects, which must be studied is the time 
of quark matter nucleation. It is still not clear whether such transition takes 
place immediately after the supernova explosion or much latter [13]. If it occurs 
(if at all) almost simultaneously with the explosion, then what is the effect of 
strong magnetic field present in the neutron star on the (first or second order) 
quark-hadron phase transition? It is also interesting to study the effect of strong 
magnetic field on the equation of state of quark matter. 

Now from the observed features in the spectra of pulsating accreting neutron 
star in binary system, the strength of surface magnetic field of a neutron star is 
found to be ~ 10 12 G. In the interior of a newly born neutron star, it propably 
reaches ~ 10 18 G [14]. It is therefore advisable to study the effect of strong magnetic 
field on SQM, which we expect to present at the core of a compact neutron star. As 
is well known, the energy of a charged particle changes significantly in the quantum 
limit if the magnetic field strength is equal to or greater than some critical value 
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Bm = mfc 3 / (qih) in Gauss; where rrii and qi are respectively the mass and charge 
(absolute value) of the particle (e.g. qi = 2e/3 for -u-quark, e/3 for d and s quarks 
and e for electron, here e =| e | is the absolute value of electronic charge), Ti and c 
are respectively the reduced Planck constant and velocity of light, both of which 
along with the Boltzmann constant ks are taken to be unity in our choice of units. 
For an electron of mass 0.5MeV, the strength of this critical field as mentioned 
above is B^^ ~ 4.4 x 10 13 G, whereas for a light quark of current mass 5MeV, 
this particular value becomes = 10 2 x B„' e ', on the other hand for s-quark of 
current mass 150MeV, it is ~ 10 20 G, which is too high to realize at the core of a 
neutron star. Therefore the quantum mechanical effect of neutron star magnetic 
field on s-quark may be ignored. The critical magnetic field as defined above is 
the typical strength at which the cyclotron lines begin to occur, and in this limit 
the cyclotron quantum is of the order of or greater than the corresponding rest 
energy. This is also equivalent to the requirement that the de-Broglie wavelength 
is of the order of or greater than the Larmor radius of the particle in the magnetic 
field. 

In recent years a lot of interesting work have been done on the properties 
of dense astrophysical and cosmological matter in presence of strong magnetic 
field [15-17]. How the primordial magnetic field affects the expansion rate of the 
Universe and synthesis of light elements in the early Universe micro-second after 
Big Bang has recently been done by Schramm et al [18] (see also ref. [19]). In 
a series of publication we have also reported the work done by us on the effect 
of strong magnetic field on SQM using both MIT bag model as well as D3QM 
model of confinement [20-23] . We have also done some work on the magnetostatics 
of superconducting quark stars and also on the stability and gross properties of 
strange stars in presence of strong magnetic field [24] . 

4 



The aim of the present paper is to investigate the stability and some gross 
properties of quark matter in presence of strong magnetic field. In section 2 we 
have studied the thermodynamic properties of bulk SQM in presence of strong 
magnetic field. In section 3 the thermal nucleation of stable quark matter bubbles 
in compact metastable neutron matter at the core of a neutron star in presence 
of strong magnetic field has been investigated. In section 4 the properties of 
co-existing bulk phases in presence of strong magnetic field has been studied. 
Assuming the possibility of a metal-insulator type of second order quark-hadron 
phase transition at the neutron star core at zero temperature, we have investigated 
the effect of strong magnetic field on such transition in section 5. In section 6, the 
paramagnetic properties of SQM has been discussed in detail. Whereas, section 7 
contains conclusions of this work. 

2. Bulk SQM in Strong Magnetic Fields 

To study the deconfining transition at the core and the stability of SQM in 
presence of strong magnetic field, we have considered the conventional MIT bag 
model. The limitation of this model has already been discussed in the introduction. 
For the sake of simplicity we are assuming that quarks are moving freely within 
the system and as usual the current masses of both u and rf-quarks are extremely 
low (in our actual calculation we have taken current mass for both of them to be 
5MeV, whereas for s-quark, the current mass is taken to be 150MeV). 

For a constant magnetic field along the z-axis (B m = B m ^ = B m = con- 
stant) , the single particle energy eigen value is given by [25] 

41* = [k 2 + + q t B m (2n + s + 1)] 1/2 (2.1) 

where n = 0,1, 2, being the principal quantum numbers for allowed Landau 
levels, s = ±1 refers to spin up (+) or down (— ) states and k is the component of 
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particle momentum along the direction of external magnetic field. Setting 2u = 
2n + s + 1 , where v = 0,1, 2, we can rewrite the single particle energy eigen 
value in the following form 

c « = [k 2 + m 2 l + gi B m 2u] 1/2 (2.2) 

Now it is very easy to show that v = state is singly degenerate while all other 
states with v ^ are doubly degenerate. 

The general expression for thermodynamic potential of the system is given by 

Q = -TlnZ 

= Pi,v(T, fn)V + Qi, s (T, fii)S + Qi,c(T, ^)C] + BV (2.3) 

i 

where the sum is over u, d, s-quarks and electron (e). 

The first term on the second line of the above equation is the volume contri- 
bution, whose explicit form is given by [26] 

fi i>v (T,^) = -7^3 J d 3 k\n(l + eM/3(^-^))) (2-4) 

where gi is the degeneracy of the i-th species (= 6 for quarks and 2 for electron) 
and V is the volume occupied by the system. 

The second term comes from the surface contribution and is given by [27] 

~ / x grT f d 3 k \ 2 

n^^i^y^! __»»-. (_) 

ln(l + exp(/?(Aii-ei))) = <7 (2.5) 

and S is the area of the surface enclosing the volume V and a is the surface energy 
per unit area or the surface tension of quark matter. 

The last term corresponds to curvature correction, which becomes zero for 
a plane interface separating two phases. The explicit form for curvature term is 
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given by [28] 



a,c(T, fJLi) 
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(2.6) 



where C is the length of the line element drawn on the surface S. 

At finite T ~ fa one has to consider the presence of anti-quarks and positrons 
in the system. We assume that for anti-particles fa = —fa, i.e. they are also in 
chemical equilibrium with respect to annihilation processes. Now we would like 
to see the necessary changes for these expressions which has to be incorporated, 
if a strong external magnetic field is present in the system. With these modified 
form of thermodynamic potentials, we shall investigate the magnetism arising from 
quantization of orbital motion of charged particles in presence of strong magnetic 
field, known as Landau Diamagnetism. It is well known that in presence of a 
magnetic field along z-axis, the path of the charged particle will be a regular helix 
whose axis lies along the z-axis and whose projection on x — y plane is a circle. 
If the magnetic field is uniform, both the linear velocity along the field direction 
and the angular velocity in the x — y plane will be constant, the latter arises from 
the constant Lorentz force experienced by the particle. Quantum mechanically 
the energy associated with the circular motion in the x — y plane is quantized 
in units of 2qiB m . the energy associated with the linear motion along z-axis is 
also quantized; but in view of the smallness of the energy intervals, they may be 
taken as continuous variables. We thus have eqn.(2.1) or (2.2) as single particle 
energy eigen value. Now these magnetized energy levels are degenerate because 
they result from an almost continuous set of zero field levels. All these levels for 
which the values of the quantity k 2 x + ky lies between 2qiB m u and 2qiB m (u + 1) 
now coalesce together into a single level characterized by the quantum number v. 
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The number of these levels is given by 



S 1 f dk x dk y = ^^ (2.7) 



(2tt) 2 J J x y 2tv 

here S is the area of the orbit in the x — y plane. This expression is independent 
of v. Then in the integral of the form f d 3 kf(k), we can replace f f dk x dk y by the 
expression given above, whereas the limit of k z , which is a continuous variable, 
ranges from — oo to +00. Then we can rewrite the eqns.(2.4) to (2.6) in presence 
of strong magnetic field in the following form 

r> 00 />oo 

tli,v = -T^f^ £ / dk * ln f 1 + e MP(»i ~ 4 v) )j) (2-8), 



^LS = T 



qi9iB r 
16tt 
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and 



n 
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dk 7 



ln (l + exp^-e^))) 



1 - 
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TV 
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rrii 



where k^ = 2uqiB m 

Since the surface and curvature terms play significant roles during quark bub- 
ble nucleation in dense neutron matter, and are not at all important for a bulk 
quark matter system, we shall study their dependence on strong magnetic field in 
the next section, when we shall discuss thermal nucleation of quark droplets in 
presence of strong magnetic field. In this section we shall consider only the bulk 
term. 



(2.10) 



We shall use eqn.(2.4) for s and s and eqn.(2.8) for u, d and e and their anti 
particles. In eqn.(2.4), = (k 2 + m 2 ) 1 / 2 , whereas in eqn.(2.8), it is given by 
eqn.(2.2). 

Then from the well known thermodynamic relations [26], the expression for 
kinetic pressure of the system is given by 

P = -J2^,v (2.11), 

i 

the number density of the ith species is given by 

--(^X (2 - i2) 

and finally the expression for the corresponding free energy density is given by 

Ui = Qi, v + mm - T (2.13) 

The last term in eqn. (2.13) comes from the non-zero entropy of the system, which 
is zero for T = 0. Then the total energy density of the confined SQM is given by 

U = J2 U i + B (2-14) 

i 

To evaluate all these thermodynamic quantities as defined above for a given 
temperature and bag parameter, we need chemical potentials for all the con- 
stituents present in the system. Assuming the condition of /3-equilibrium, we 
have 

Hd = Hs= H (say) (2.15) 

and 

li u = [i-li e (2.16) 
Whereas the charge neutrality condition gives 

2n u — rid — n s — 3n e = (2-17) 
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Again the baryon number density of the system is given by 

riB = \( n u + n d + n a ) (2.18) 

which is considered as a constant parameter. 

We have solved these equations numerically to obtain chemical potentials for 
all the specieses for various values of B m (including zero) and then evaluate the 
thermodynamic quantities of the SQM system. 

Since for any non-zero temperature T » 0, analytical expressions for the 
thermodynamic quantities can not be obtained, we shall give zero temperature 
and for the sake of illustration low temperature analytical expressions. Whereas 
in actual calculation for T ^ 0, we shall use numerical solutions of the general 
expressions. 

Now for T = 0, we have the number density of the component % {% = u,d or 

e) 



m = n i 



max J 



c 1/2 + 2 y.^--) 1 ' 2 



(2.19) 



where 



n i = ^J2^9i (2-20) 

Q = -< (2 21) 

2(/i B m 

and [^mij is the greatest integer not exceeding Q. The v limit becomes finite for 
T — > limit only, when the maximum available energy of a particle ~ its Fermi 
energy. On the other hand, for s-quark, the expression for number density is given 
by the usual expression 

n s = 5 (2.22) 

In fig.(l) we have shown the variation of number densities for u, d, s-quarks 
and electron with the magnetic field intensity B m . At low enough magnetic field 
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strength, all the number densities remain almost constant and are equal to their 
zero field values. For relatively higher field strengths > 10 3 x B$ G, electron 
density increases with B m . The s-quark density decreases with B m at relatively 
higher values than that of electron and becomes negligibly small for extremely 
large field strength > 10 5 x B^^G. On the other hand, u quark density be- 
comes almost two times the zero field value at this extremely large field strength. 
Whereas the <i-quark density remains almost constant within this range of mag- 
netic field. But all the number densities show oscillating behavior with the increase 
of magnetic field intensity. The number densities show oscillating nature, as the 
consecutive Landau levels for different constituents pass through the correspond- 
ing Fermi levels. The electron density increases to ~ 10 4 times zero field value at 
extremely large magnetic field. One can very easily check that these oscillating 
and increasing and / or decreasing natures of number densities are consistent with 
the charge neutrality and constancy of baryon number density conditions. The 
physical meaning of change in s-quark density and electronic density is the shift 
of /3-equilibrium point in presence of strong magnetic field. As a cross check one 
should try to see how do the rates of weak processes in which s-quark and / or 
electron are created or annihilated, change in presence of strong magnetic field (a 
detail analysis on the effect of strong magnetic field on weak and electromagnetic 
processes in SQM will be reported somewhere else [29]) 

The expression for kinetic pressure and free energy density for the ith species 
are given by (i = u,d and e) 
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where = (mf + 2uq i B Tn ) 1 ^ 2 . Whereas for s-quark, they are given by the usual 
expressions 



P, = - mlf /2 - Sm^.tji 2 , - ml) 1 ' 2 

+ 3m1lni"- + ^-"'-> 1/2 }] (2.25) 
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^ s + (^ 2 -m 2 ) 1 /2 



(2.26) 



These two expressions are without the bag pressure or vacuum energy density. 
Here we have also not taken magnetic pressure and magnetic energy density into 
account. They are given by the usual expressions, ±5^/8^, where + sign is for 
magnetic energy density and — sign is for magnetic pressure. Now for B m = 10 15 G, 
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the magnetic energy density ei 5 = 32.283MeV 4 , which is too small compared 
with the free energy density ~ 10 6 MeV 4 , whereas for B m ~ 10 18 G, it is given 
by ei8 = 10 6 x 615. The same conclusion is also valid for magnetic pressure of 
the system. Therefore, magnetic energy density as well as the magnetic pressure 
will be comparable with the particle free energy density and the kinetic pressure 
respectively if and only if the magnetic field intensity is too high. For B m > 10 20 G, 
the magnetic energy density as well as the magnetic pressure play the main roles. 
This magnetic field strength is however, too high to realize at the core of neutron 
star or any other stellar objects of our interest. 

To study the stability of bulk SQM in presence of strong magnetic field we 
shall use the stability condition J2i Pi — B = 0, where % = u,d,s and e, and solve 
the eqns.(2.15)-(2.18) numerically for T — 0. Since Vmax is a function of /Uj, we 
have solved these equations self consistently for B m 7^ 0. In fig. (2) we have plotted 
energy per baryon of SQM against the baryon number density. The upper curve 
is for B m = 0, whereas the lower one is for B m = 10 3 x B„ ^ e '. This figure shows 
clearly that the presence of strong magnetic field makes SQM energetically more 
stable. 

We shall now write down for the sake of completeness the low temperature an- 
alytical expressions for thermodynamic parameters. The number density is given 



for i = u,d and e, and as before, for s-quark it is given by 

s ~ s + 6 (^_ m 2)l/2 I' 
The expressions for pressure for u, d and e is given by the general expression 



by 




(2.27) 




(2.29) 



13 



and for s-quark it is given by 

Ps = P^ + ^M-ml)^ (2.30) 

In all these eqns., the supescripts "0" refer to zero temperature cases. In this 
particular low temperature case the free energy density for the ith species is given 
by 

e, = -Pi + fi.n, + Ts z (2.31) 

where i = u, d, s or e. here the last term comes from non-zero entropy of the 
system and is given by 

(CIO \ O m ax * 

for % = u, d and e, whereas for s- quark it is given by 

Tss = ~t(^\ =T^M-mlY" (2.33) 

In figs. (3) and (4) we have shown the equation of states for SQM at T = 
and 40MeV respectively. In both these figs, the dotted curves are for B m = 0, 
whereas the solid ones are for B m = 10 3 x B™" 6 '. For the sake of comparison we 
have plotted B-J equation of state [14] in fig. (3) and is indicated by the symbol bj. 
These two figs, show that in presence of strong magnetic field the equation of state 
of SQM changes significantly. In both these figs, we have taken B 1 / 4 = 0, which is 
a constant scaling factor. Fig. (5) shows the variation of energy per baryon of stable 
SQM with the magnetic field strength for n B = 2.5n and T = 0. This fig. shows 
that for low B m , the energy per baryon remains almost constant and rises sharply 
as B m exceeds 10 4 x B^^, and the matter becomes energetically unstable. In 
fig. (6) we have plotted the variation of kinetic pressure of SQM with the magnetic 
field intensity for B 1 ^ 4 = and for the same value of ns as before. This fig. shows 
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that for relatively low values of magnetic field strength the kinetic pressure initially 
remains almost independent of magnetic field strength and then decreases sharply 
as B m exceeds 10 4 x Bm' 6 ' and finally becomes negative for extremely large values 
for B m . Therefore for extremely large magnetic field strength, the system becomes 
energetically as well as mechanically unstable (these two conclusions are valid if 
the magnetic field is is confined within the quark matter region only, which is 
possibly not true). 

3. Nucleation of Quark Bubble in Compact Neutron matter 

In this section we shall investigate the first order quark-hadron phase transi- 
tion initiated by the nucleation of droplets of quark matter in presence of strong 
magnetic field. In particular we shall study the effect of strong magnetic field on 
the nucleation of quark droplets in dense neutron matter. Since the surface and 
curvature energies of the quark bubble play crucial role in droplet nucleation we 
would like to see how these two quantities are modified by the presence of strong 
magnetic field. 

Now the nucleation rate of stable quark matter bubble in metastable neutron 
matter per unit volume due to thermal fluctuation is given by [30] 

I = I exp(-W m /T) ^T 4 e*p(-W m /T) (3.1) 

where W m is the minimum thermodynamic work to be done to create a critical 
quark droplet and is given by [31] 

Wm = l n JKPY [2 + 2(1 + 6)3/2 + 3b] (3 - 2) 

where a = a q + a n the surface tension and AP = P q — P n is the pressure difference. 
Here AP = P q — P n is a positive quantity and both a q and a n are also positive, 
b = 27(AP)/a 2 , and 7 = 7 g — 7 n , stands for curvature energy density. Here phase 
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"q" being a droplet of quark matter and phase "n" stands for the metastable 
neutron matter. In eqn.(3.1), the pre-exponential factor is chosen to be T 4 , where 
T is the temperature of the metastable medium (for an exact expression see ref. 
[32]). In eqn.(3.2) we have assumed that the two phases are in chemical equilibrium 



the strong interaction time scale, the creation of strange quarks through weak 
processes within the quark droplets may be ignored. Therefore the constituents 
of quark droplets are only u and d quarks. Since the temperature (~ 4 — 5 MeV) 
<< quark chemical potential (~ 300 MeV), the presence of anti-quarks can also 
be ignored. On the other hand, if we assume the presence of hyperons at the core 
of neutron star (which is believed to be true), then immediately after the phase 
transition to quark matter, s quarks will also be present in the quark bubble. 
However, the quantum mechanical effect of strong magnetic on s quark part of 
strangelet can be ignored as long as the magnetic field strength is < 10 20 G, which 
has already been discussed. Therefore in the surface energy term of strangelet 
a finite contribution will come from s quark part and is given by eqn.(2.5). On 
the other hand the surface tension or the surface energy per unit area for u and 
d quarks are given by eqn.(2.9). It is obvious from this eqn. that the surface 
energy diverges logarithmically in the infra red limit (i.e., k z — > 0) for v = 0, i.e. 
for the ground state Landau level. To show this more explicitly, we shall make 
zero temperature approximation. In this case the upper limit of the u-sum can be 
obtained from the condition 



(Afj, = 0). Since the bubble nucleation time Tbubbie 



1Q -23 



sec 



k 2 z = fJ% ~ ml - 2vqiB m > 
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This gives 
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WiSax] (say) 
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which is the largest integer not exceeding (//| — m?)/(2qiB m ). Therefore as has 
been mentioned earlier the upper limit [^L] can not be same for u and d quarks. 
Now to visualize the effect of magnetic field on the nucleation rate of quark droplet, 
we shall first rewrite eqn.(2.9) in the form 



TB m 
16n 



1^ / 77% ;2 



In 



1 + exp - 



O) 



G (3.5) 



where 



G = 







1 - - tan" 1 
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for MIT bag model and is equal to unity for D3QM model. Now we shall evaluate 
integral (3.5) by parts for G = 1 and take T — > limit, then 



B 



i=u,c 



k z dk z 



„=o Jo yjh* + m? + k\ 



\n(k z + Jk1 + k 2 Mi) ) 



v=0 



(3.6) 



It is easy to check that eqn.(3.6) is also a part of a, given by eqn.(3.5) for which 
G^l. Since \ii > \Jk^_ ^ + mf , therefore, for v = the second term of eqn.(3.6) 
becomes +oo. The other part of integral (3.6) has been evaluated numerically and 
found to be a finite number. Similarly the integral (3.5) with the second part of 
G (which is 7^ 1) has also been obtained numerically and is found to be finite. 
Therefore the diverging property of a(G = 1) for v = is also true for the whole 
surface energy / area, given by eqn.(3.5). The term v = corresponds to the 
lowest Landau level. If the magnetic field is too strong only this level will be 
populated, which is also obvious from eqn.(3.4). Although we have assumed here 
T — > 0, this important conclusion is equally valid for any finite T. Since both P n 
and P q are finite even in presence of strong magnetic field and T is also finite, we 
have from eqn.(3.1) 7 = 0, which means the thermal nucleation rate of droplet 
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formation becomes zero, i.e., there can not be a single quark droplet formation 
from metastable neutron matter due to thermal fluctuation at the core, if the 
magnetic field is of the order of or greater than the corresponding critical value. 
Therefore the formation of quark droplet at the magnetized neutron star core will 
be controlled by some other mechanism, namely, the quantum effects, triggering 
by strangelet capture etc. 

Following ref. [28] , we shall now investigate the effect of strong magnetic field 
on the quark curvature term which also plays an equally important role as surface 
tension during thermal nucleation. As before, we can rewrite eqn.(2.10) in the 
form 

7 = 7^2 9iQi Yl / k 2 T P ln 0- + exp(-/3(e, - &))) G (3.7) 

i=u,d u=Q J ° z -L,(«) 

where 

G=l--— - -tan" 1 — 
Zrrii 2 rrii 

for MIT bag model and equal to unity for D3QM model. 

Now we shall evaluate eqn.(3.7) term by term. Consider 

r°° dk 

h = / u2 ■ A ln C 1 + ex P(-/3(^ "/*))) ( 3 - 8 ) 

Jo K z -t- 

Obviously, this integral diverges in the infra red limit for v = 0. Here the diver- 
gence is l/k z type. To see it more explicitly, let us integrate eqn.(3.8) by parts, 
then we have 

Il = TkI^l ^ 1{ kZ^ ] {kl + k{ Ai) + m?)V2 exp(/3(e l - + 1 (3 ' 9) 

which diverges for v = 0, but the divergence is not logarithmic (ii ~ 1/u as 
v -> 0). 

18 



Next consider the second term 
^7r f°° elk 

Which has logarithmic divergence in the infra red limit for v = 0. Integrating by 
parts, we have 

h = ^-[{ln(/c ±)(i) )ln(l + exp(-/?(e i -/x i )))} 

"^1 ln(jfcz + (Jfe * + «« )1/2) (jfe2 + ki^l mty/i 
1 

_exp(/3(e; - fj,i)) + 1 

For z/ = 0, the first term diverges logarithmically, but unlike the second term of 
eqn.(3.6), it becomes — oo. The second term of eqn.(3.11) and the last term of 
eqn.(3.7) have been checked numerically and they remain finite for all values of 
v in the infra red as well as ultra violet limits. Now the divergences of first two 
terms of this eqn. can not cancel each other. The first divergence is much faster as 
v — > than the second one, therefore the overall divergence of curvature term (3.7) 
remains positive as v — > 0. Now from eqn. (3.2) it is obvious that if some how the 
quantity b becomes finite (even zero), still the effect of a makes W m infinitely large. 
Therefore, in general, both the surface tension as well as curvature terms of a quark 
matter bubble diverge in presence of strong magnetic field. Physically it means, 
that an infinite amount of work has to be done to create a critical quark bubble. 
As a consequence, thermal nucleation of such bubble in metastable neutron matter 
will be completely forbidden. It can only occur by some non-thermal mechanisms. 
Therefore, in future, if we get some direct or indirect experimental evidences for 
quark core and if the above conclusion is assumed to be correct, then either the 
magnetic field of such objects will be extremely low (< Bm ) from its time of birth, 
if quark matter is assumed to be produced immediately after supernova explosion, 
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(3.11) 



or it must have been produced much latter, when the neutron star magnetic field 

(c) 

has decayed to a value B m < Bm or the quark droplets are produced through 
some non-thermal means. Since all these divergences are appearing in the infra 
red limit for v = 0, one can in principle remove them by introducing a lower cutoff 
for k z . Ofcourse, in that case one has to think the physical meaning of this lower 
cutoff of longitudinal momentum. 

4. Co-existing Bulk Phases 

From the discussion of previous section, we have seen that it is impossible 
to have a mixed phase of dense neutron matter and droplets of quark matter in 
presence of strong magnetic field of strength greater than the corresponding critical 
value. In this section we would like to see the effect of strong magnetic field on 
the coexisting bulk phases. 

Assuming for the sake of simplicity that only neutrons present in the hadronic 
phase and obey the free hadronic gas equation of state; then the conditions to be 
satisfied by the co-existing phases are 

A*n = A*« + 2/Lid (4.1) 

Pn = P q (4.2) 

and 

T n =T q = T c (4.3) 

We also assume that the condition of /3-equilibrium is satisfied in bulk quark matter 
phase. 

In fig. (7) we have shown the phase diagram (n B vs T diagram) for such a 
mixed phase. The dotted curve is for B m = 0, whereas solid curve is for B m = 
10 3 x Bt )(e) . This figure shows that the phase diagram changes significantly in 
presence of strong magnetic field. 
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This observation is specially important during quark-hadron phase transition 
in the early Universe (where tib ~ 0) if it would have taken place in presence of 
strong magnetic field (!!) of our interest (for a detail discussion, see the in press 
paper of ref. [23] , where the effect of strong magnetic field on the baryon number 
inhomogeneity during cosmic QCD phase transition has also been discussed). 

5 Metal-Insulater Type of Transition 

Since the order of quark-hadron trasition is still not known exactly, in this 
section we shall consider a metal insulator type of second order quark-hadron 
phase transition at zero temperature in presence of strong magnetic field. Now as 
the density of neutron matter at the core of neutron star increases, the separation 
between two neutrons decreases and at some critical value, when the separation 
becomes ~ 0.4fm, the hard core radius, neutrons virtually loose their individuality 
and quarks can percolate out, give rise to a color metal. This transition is more 
or less identical with the metal insulator transition observed in condensed matter 
physics, which takes place at high pressure. Here neutrons, which are colorless 
objects play the role of neutral atoms of the insulator placed at the lattice points 
and orbital bound electrons may be compared with the confined quarks. 

Unlike a first order transition, where the two phases are in equilibrium at 
the critical point, here they are indistinguishable. The conditions which must be 
satisfied in this case at the critical point are 



n 



(5.1) 



= 2\± d + fi. 



u 



(5.2) 



and of course 



P = P 

1 n 1 q 



(5.3) 
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Here we have assumed the presence of u and d quarks only, the numerical 
solution of the eqns.(5.1-5.3) give critical baryon number density n c = 1.49no and 
S 1 / 4 = 146.52MeV for B m = 0. But in presence of strong magnetic field, both 
the critical baryon number density and the corresponding bag pressure depend on 
the strength of magnetic field. In fig. (8) we have shown the variation of n c with 
the magnetic field intensity B m and in fig. (9) the same type of variation is shown 
for B 1 / 4 . These figs, show that for low magnetic field strength, both the critical 
density and bag pressure are almost independent of B m and very close to their 
zero magnetic field values. But for large magnetic field strength > 10 4 x Bm^ e 'G, 
both these quantities increase sharply. These figs, also show that for extremely 
large magnetic field strength, the critical density becomes too high to achieve at 
the core of neutron star. Which implies that for an extremely strong magnetic field 
at the neutron star core, both the first order transition through the nucleation of 
quark bubbles and a second order metal insulator type of quark hadron transition 
are impossible. 

6. Pauli Paramagnetism of Quark Matter 

As has been mentioned in the introduction, that the magnetism of quark gas 
at T = consists of two parts: a paramagnetic part due to the intrinsic spin 
magnetic moment, known as Pauli paramagnetism and a diamagnetic part called 
Landau diamagnetism, arises from the quantization of the orbital motion of the 
quarks in the strong magnetic field. The second part has already been discussed. 
In this section we shall consider the Pauli paramagnetism of quark matter. 

It is known that the quarks have non-zero magnetic dipole moment, given 
by T] u = 1.852rjN, i]d = —0.972^ and r] s = — 0.581^, where tjn = e/2m p = 
3.152 x 10~ 14 MeV/Tesla, the nuclear magneton. Whereas for electron, it is given 
by 7] e = 1.0011677s, where r] B = e/2m e = 5.788 x Kr n MeV/Tesla, the Bohr 
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magneton. 

Therefore in presence of a magnetic field of strength B m , the magnetic poten- 
tial energy of the ith component is —r\iB m . Where % stands for u, d or s-quarks 
or electron. Then the corresponding number density is given by [33] 



n 4 



9% 
2tt 2 



k 2 dk+ / k 2 dk 

e=m,i J e=m i +27j i B m 



which gives after some simple algebraic manipulation 



n. 



_ 2r)iB m (2i]iB m + 2mj) (n 3/2 _ ^ 
6yr 2 K^i V 



where 



X 



e 2 (B m )-m 2 



(6.1) 



(6.2) 



(6.3) 



2rj i B m (2r]iB rn + 2m*) 
and is the single particle energy. From eqn.(6.3) one can write down 

e 2 (B m ) =m 2 + x z [R 3/2 ]- 2 / 3 (el - m\) (6.4) 

where e^o is the single particle energy in absence of magnetic field, given by 



Ci,o = m 2 + (37T 2 ?!,) 2 / 3 



(6.5) 



and 



R 3/2 ( Xi ) = ^(2x s / 2 -l) 



<:; 2(-'- = -I--'' - 1 ; ( 6 - 6 ) 

Now the total particle number density for the ith component can be decompossed 
into two parts; with the magnetic moment along the direction of magnetic field 
and opposite to it, then we have 



n i = n? ) (B m )+n\ i \B m ) 



(6.7) 



where 



n^(B m ) = ^ 2 [R 3/2 (x l )]^ 
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(6.8) 



the number density for the ith component with the direction of magnetic dipole 
moment along the direction of external magnetic field and 

n^\Bm) = y^Wo^r'a - *f' 2 ) ( 6 - 9 ) 

the same number density as above, but the direction of magnetic dipole moment 
is now opposite to the direction of external field. In absence of magnetic field, or 
for B m — > 0, the quantity Xi = oo, which gives 

e i = e ifi , n { P = n { t l) (6.10) 

The second relation implies equal occupation probability for up and down spin 
states. On the other hand for Xi = 1, we have 



n { p = m, and nj i} = (6.11), 



which is the complete saturation of the spin states along the direction of magnetic 
field. In this condition 

ef = mf + 2 2 / 3 (e2 -mf) (6.12) 

The threshold magnetic field strength for complete saturation of ith component 
can be obtained by putting Xi = 1 in eqn.(6.3), which gives 



B (th)(i) _ 



22/ 3 ^y /3 +m 2 



1/2 



- mi 

(6.13) 



2Vi 

For rii = 2.5ns, the threshold values for it, d and s quarks are 3.4 x 10 19 G, 6.46 x 
10 19 G and 7.59 x 10 19 G, respectively, which are of course not too low. 
The energy density for such a Pauli paramagnetic system is given by 



U = Y, u i ( 6 - 14 ) 
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where 



Ui(B m ) = 



9% 
2tt 2 



ek 2 dk + / ek 2 dk 

e=m; J e=m i +2r]iB rn 



(6.15) 



which after some straight forward algebraic manipulation becomes 



Ui(B m ) = -|4 [2(e 2 - m 2 ) 1 / 2 (2e 2 - m 2 )e, - 2m 4 In 



2 W 2 f9, 2 



16tt 2 

- (yt(l/i + 2mi))(m? + 5m 2 ?/; + 6m;?/ 2 + 2yf) 



e, + (e 2 -m 2 )V2 



+ mfln(l + K + fe(2 ^ + ! " ))1/2 )] 



where ^ = 2^S m 

The kinetic pressure is given by 



(6.16) 



(6.17) 



where 



Pi(B m ) = 



9% 
6tt 2 



e(B m ) 



£; 4 d£; 



+ 



fc 4 d£; 



(fc 2 + m 2 )V2 y e=mi+27?i(fc 2 + m 2 )1/2 
after some algbraic simplification it reduces to 

Pi(B m ) = ^ [2(e 2 - m 2 )V 2 (2e 2 - m 2 )e t + 6mf In ( ^^f^ ) 



(6.18) 



- + 2m i ))(-3mf + m 2 ?/; + 6m,?/ 2 + 2y?) 

3 m ?lnfl+ K + fa(2 "" + y ' ))1/2 )l 



(6.19) 



In the field free case (B m = 0), y^ = 0, then both the expressions for energy density 
and pressure reduce to the corresponding zero field expressions, with e; replace by 
e\°\ On the other hand for the saturation condition, X{ = 1 and yi = — m t , 
then we have 



C^a* = 1) = ^ [2(e 2 - m 2 ) 1 / 2 (2e 2 - m 2 )e, - 2mj In 



4i„ / e i + (ei-™>i 



2U/2 



(6.20) 
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where in this expression ej is given by eqn.(6.12). Similarly the expression for 
pressure is given by 



9i 



Now the bulk interaction energy per unit volume is given by 



(0 

in 



(6.21) 



(6.22) 



where 



(6.23) 



Substituting the values of and and after some simple algebraic manipu- 



lation, we have 



U-n = niTJiBr 



3/2 
1 X i 



Rz/2{Xi) 



(6.24) 



which becomes zero for B m = and at the saturation limit, 



U$(Xi = 1) = UiTjiBr 



(6.25) 



Finally, the chemical potential for the species i in presence of a magnetic field 
B m is given by 

m = ei- m B m (6.26) 

Substituting the values for ci and rjiB m , we have after some simple algebra 

1/2 l 



Hi = [mf +Xi[R 3/2 ] 2/3 (4 >0) -m, 2 )} " - - 



ef -m?(l -X?) 1 1/2 



X 7 ; 



- m, 
(6'.27) 

which becomes e^o for y$ = 0, or equivalently for £? m = or — > oo. On the 
other hand in the saturation limit, we have 



(jLifai = 1) = - [et(xi = 1) + m i 



(6.28) 
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7. Conclusions and Discussions 

From the theoretical investigations of SQM in presence of strong magnetic 
field we can make the following conclusions: 

(i) In presence of a strong magnetic field of strength of the order of greater than 
some critical value as discussed in the main text, the SQM becomes energetically 
more stable as long as B m < 10 20 G. Above this value since the magnetic energy 
density and also the magnetic pressure play the main roles the system becomes 
energetically as well as mechanically unstable (pressure of the system becomes 
negative). 

(ii) Equation of state of magnetized SQM changes significantly. The (3 equi- 
librium point has also been shifted. 

(iii) There can not be any thermal nucleation of quark droplets in a metastable 
state of compact neutron matter. 

(iv) If bulk QGP and hadronic matter can co-exist in presence of strong 
magnetic field, then a significant change in phase diagram has also been observed. 

(v) If the quark-hadron phase transition at low temperature and high density 
is like a metal insulator type of second order transition. Then in presence of 
extremely strong magnetic field as mentioned above, the critical density at which 
such transition may occur is found to be too large, which is almost impossible to 
achieve at the core of a neutron star. 
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Figure Captions 

1. The variation of u, d, s-quarks and electron number density with the magnetic 
field intensity B m , ub = 2.5no- 

2. Plot of energy per baryon against the baryon number density for B m = 
(upper curve) and B m = 10 3 x B^^ (lower curve). 

3. Pressure vs energy per baryon curves (equation of states) for T = 0. The 
dotted curve is for B m = and solid one is for B m = 10 3 xBl c)(e) . The curve 
indicated by bj is for B-J equation of state. 

4. Pressure vs energy per baryon curves (equation of states) for T = 40MeV. 
The dotted curve is for B m = and solid one is for B m = 10 3 x B™" 6 '. 

5. variation of energy per baryon for stable SQM with the magnetic field strength 
for ub = 2.5no and T = 0. 

6. variation of pressure with the magnetic field strength for n B = 2.5n , T = 
and B 1 / 4 = 0. 

7. Phase diagram for a bulk co-existing QGP and hadronic phases for B m = 
(dotted curve) and B m = 10 3 x B™' 6 ' (solid curve). 

8. variation of n c with the magnetic field intensity B m . 

9. variation of S 1 / 4 with the magnetic field intensity B m . 
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